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Abstract 



We report on the measurement of the direct emission (DE) and interference (INT) 
terms of the tt^tt^j decay by the NA48/2 experiment at the CERN SPS. 

From the data collected during 2003 and 2004 about 600k such decay candidates 
have been selected. The relative amounts of DE and INT with respect to the inter- 
nal bremsstrahlung (IB) contribution have been measured in the range < T* < 
80 MeV: 

PracDE(0 < T* < 80 MeV) = (3.32 ± 0.15stat ± 0.14,,,,) x 10-^, 
FraciNT(0 < T; < 80 MeV) = (-2.35 ± 0.35,tat ± 0.39,^,) x lO^^ 

where T* is the kinetic energy of the charged pion in the kaon rest frame. This is 
the first observation of an interference term in — )• tt^tt^'-^ decays. 
In addition, a limit on the CP violating asymmetry in the and K~ branching 
ratios for this channel has been determined to be less than 1.5 xlO^'^ at 90% 
confidence level. 

1 Introduction 

The decay channel — ir'^TT^'-f is one of the most interesting channels for studying 
the low energy structure of QCD. Radiative nonleptonic kaon decays, such as Kl,s 
7r+7r^7 and — t- tt^tt^ are dominated by long distance contributions, but a small short 
distance contribution is also present and can be studied. 

The total amplitude of the K"^ — )■ ir'^'n'^'-f decay is the sum of two terms: the inner 
bremsstrahlung (IB) associated with the — )■ vr^vr^ decay with a photon emitted from 
the outgoing charged pion, and the direct emission (DE) in which the photon is emitted 
at the weak vertex. Using the Low theorem [1] the branching ratio of the IB component 
can be predicted from that of the — > vr^vr'' channel, using QED corrections [21 [3]. 
As the K"^ — )■ tt^tt" decay is suppressed by the A J = 1/2 rule, the IB component of 
— )■ 7r^7r°7 will be similarly suppressed, resulting in a relative enhancement of the DE 
contribution. 

The DE term has been extensively studied in the framework of Chiral Perturbation 
Theory (ChPT) |H El El [3 [H |9] . Direct photon emission can occur through both electric 
and magnetic dipole transitions. The electric dipole transition can interfere with the 
IB amplitude giving rise to an interference term (INT), which can have CP violating 
contributions. In ChPT, DE arises only at order O(p^) and cannot be evaluated in a 
model independent way. The magnetic part is the sum of two anomalous amplitudes: 
one reducible, that can be calculated using the Wess-Zumino-Witten functional fiO\ ITT], 
and one direct amplitude, whose size is not model independent but is expected to be 
small. The electric transition amplitude depends on unknown constants that cannot be 
determined by symmetry requirements alone. An experimental measurement of both DE 
and INT terms allows the determination of both the electric and magnetic contributions. 

The properties of the — > Tc^Tc^'-f decay can be conveniently described using the 
T*, W variables, where T* is the kinetic energy of the charged pion in the kaon rest frame 
and ly is a Lorentz invariant variable given by [21 E] : 

here Pk-, Pn, -P7 are the 4- momenta of the kaon, the charged pion and the radiative photon. 
Values of W can vary within the range < 1^ < 1. 
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Using these variables, the differential rate for the K"^ — ti'^tx^'-) process can be 
written as [2[ |3] : 

mt,nUXl + Xl,)W'] , (2) 

where Qrp*Q^ is the differential rate for the IB component, is the CP violating phase, 
Sj are the strong pion-pion re-scattering phases for a final state of isospin I and orbital 
momentum / of the tttt system, and Xe, Xm are normalized electric and magnetic ampli- 
tudes respectively, which are W-independent quantities. Recently, the presence of a form 
factor in the pole part of the magnetic amplitude has also been suggested [12] . 

The DE term is proportional to and the INT term is proportional to W^. 
This different W dependence allows the extraction of the different decay components. In 
particular the ratio of the data W distribution with respect to a simulation of the Inner 
Bremsstrahlung (MC(IB)) component can be parameterized as: 

Data/MC{IB) = c(l + (a ± e)W^ + hW^). (3) 

The c parameter represents an overall normalization factor, and the a, h and e coefficients 
are related to the branching ratios of Direct Emission and Interference respectively by: 



RP {' dT* f W^^^^dW T 



_ hi "^^"'^ = hlEE. (4) 

dT'dW^^ _ 1 1 NT /„N 

The coefficient a parameterize the CP conserving part of INT term while the 
coefficient e is the CP violating one. As it is expected that a ^ e, the CP violat- 
ing contribution will be neglected in the study of the Interference term. In section 17.21 
the possibility of e being non zero will be investigated. Most of the previous experi- 
ments [131 HH US [El dZl UHl [19] have measured DE and INT terms in the kinematical 
region of 55 MeV < < 90 MeV, obtaining a value for the INT contribution compat- 
ible with zero. Therefore, the values quoted in these works for the DE branching ratio 
have been obtained from fits where the INT term has been set to zero. In the Particle 
Data Group (PDG) tables [20j only the latest experiments have been taken into account, 
resulting in BRde = (4.3 ± 0.7) x 10"^ 



2 NA48/2 beam line and detector 

The beamline of the NA48/2 experiment was specifically designed to measure charge 
asymmetries in — )■ tt^tttt decays [21], using secondary kaon beams produced by 
400 GeV/ c protons from the CERN SPS accelerator impinging on a Beryllium target. 
The two simultaneous oppositely charged kaon beams with central momenta of (60 ± 3) 
GeV/c are selected by a system of dipole magnets forming a so-called "achromat" with 
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Figure 1: Schematic lateral view of the NA48/2 beam line (TAXI 7, 18: monitored beam 
dump/collimators used to select the momentum of the and K~ beams; FDFD/DFDF: 
focusing set of quadrupoles, KABESl-3: beam spectrometer stations not used in this 
analysis), decay volume, and detector (DCHl-4: drift chambers, HOD: hodoscope, LKr: 
EM calorimeter, HAG: hadron calorimeter, MUV: muon veto). Note that the vertical 
scales are different in the two parts of the figure. 



null total deflection, followed by a set of focusing quadrupoles, muon sweepers and colli- 
mators (Figured]). With 7 x 10^^ protons per burst of ~ 4.5 s duration impinging on the 
target, the positive (negative) beam flux at the entrance of the decay volume is 3.8 x 10^ 
(2.6 X 10'') particles per pulse, of which 5.7% (4.9%) are Downstream of the 

second achromat, both beams follow the same path, entering a decay volume housed in a 
114 m long vacuum tank with a diameter of 1.92 m for the first 66 m, and 2.4 m for the 
rest. The beams are steered to be collinear within ~ 1 mm in the entire decay volume. 

A detailed description of the NA48 detector can be found in [22]. The charged 
decay products are measured by a magnetic spectrometer consisting of four drift chambers 
(DCH) with a dipole magnet placed in the middle. Each octagonal shaped chamber has 
4 views of sense wires, one horizontal, one vertical and two along each of two orthogonal 
45° directions. The spectrometer is located in a tank filled with helium at atmospheric 
pressure and separated from the decay volume by a thin (0.0031 radiation lengths, Xq) 
Kevlar window. Downstream of the Kevlar window the beam continues in vacuum through 
a aluminium beam pipe of 152 mm outer diameter and 1.2 mm thick, traversing the 
center of the spectrometer and all following detectors. Charged particles are magnetically 
deflected in the horizontal plane by an angle corresponding to a transverse momentum 
kick of 120 MeV/c. The spectrometer momentum resolution is (Jp/p = (1.02 © 0.044 ■ p)% 
{p in GeV/c). The spectrometer is followed by a hodoscope consisting of two planes of 
plastic scintillators segmented into horizontal and vertical strips and arranged in four 
quadrants. A quasi-homogeneous liquid Krypton calorimeter (LKr) is used to reconstruct 
7 and electron showers. It is an ionization chamber with an active volume of 7 m^ of liquid 
krypton, segmented transversally into 13248 2 cm x 2 cm projective cells by a system 
of Cu-Be ribbon electrodes, and with no longitudinal segmentation. The calorimeter is 
27 Xq deep and has an energy resolution a{E)/E = 0.032//E © 0.09/E © 0.0042 {E in 
GeV). The space resolution for single electromagnetic showers can be parametrized as ax 
= <Jy = 0A2/\/E © 0.06 cm for each transverse coordinate x, y. Due to the very good 
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determination of the shower position the calorimeter allows a very precise reconstruction 
of the 7r° mass or of the Z coordinate of the decay vertex in 7r° — )■ 77 decays. To reduce the 
data volume only information from LKr cells with a signal greater than a given threshold 
is stored. A hodoscope consisting of a plane of scintillating fibers is installed in the LKr 
calorimeter at a depth of ~ 9.5 Xq. 

NA48/2 collected data in two runs in 2003 and 2004. In order to minimize systematic 
uncertainties for the asymmetry measurements, the magnetic fields in the spectrometer 
and beam line magnets were alternated regularly. The spectrometer magnet current was 
reversed on a daily basis during 2003 and every ~ 3-4 hours in 2004. All the beam line 
magnet polarities were inverted once per week. 

For one-track events, the first level trigger (LI) requires a signal in at least one 
quadrant of the scintillator hodoscope, in coincidence with the presence of energy depo- 
sitions in the LKr geometrically consistent with more than two photons. At trigger level 
the signals from single cells are added together to form two orthogonal 4 cm wide views 
{x, y projections). 

The second level trigger (L2) is a software algorithm running on a fast on-line 
processor cluster. It receives the drift chamber information and reconstructs the momenta 
of charged particles. Assuming that the particle is a vr^ originating from the decay of a 
60 GeV/c traveling along the nominal beam axis, it evaluates the missing mass of the 
event. The requirement that the missing mass is not consistent with the 7r° mass rejects 
most of the main — > n^n^ background, reducing the rate of this trigger to ~15k 
events per burst. 

3 Event selection 

A pre-selection is performed requiring events with one charged track of momentum 
above 10 GeV/c and at least three electromagnetic clusters of energy greater than 3 GeV. 
Cluster times must be within 3 ns, and the track time within 4 ns of the mean cluster 
time. This initial sample is kept for further analysis. 

The track is considered to be a pion candidate if it lies within the spectrometer 
fiducial volume and the ratio between the energy deposited in the calorimeter and the 
momentum measured in the spectrometer {E/p) is smaller than 0.85. 

Photon candidates are defined as in-time electromagnetic clusters reconstructed 
outside a 35 cm radius disk centered on the pion impact position at the front face of 
the LKr. The number of 7 candidates must be three. To avoid too large energy sharing 
corrections, the distance between any two photon clusters must be greater than 10 cm. 
Fiducial cuts on the distance of each photon from the LKr edges and center are also 
applied to ensure full containment of electromagnetic showers. 

The closest distance of approach (CDA) between the charged pion trajectory and 
the beam axis is computed. A cut of CDA < 6 cm is required, rejecting a negligible 
amount of the signal. The position of the kaon decay vertex zch is defined at the point 
of closest approach between these two lines. 

Of the three selected photons, two of them must be associated to the decay and 
the remaining one identified as the radiated photon, also called the odd photon. Using zch 
the positions and energies of the photon clusters, the masses corresponding to the three 
possible pairings are calculated. The pairing with mass closest to the PDG 7r° mass is then 
selected as the correct 7r° pairing. For this pairing and using the PDG tt^ mass, the position 
of the decay vertex z^o is calculated and must lie in a 9000 cm long decay volume starting 
800 cm downstream from the final collimator. Only events where \zt^o — zch\ < 400 cm 
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are selected, rejecting only 1% of the signal. In addition it is required that both remaining 
photon pairings in the event satisfy the condition \zt,o — zch\ > 400 cm. This reduces the 
data sample by 20%, but also reduces the misidentification of the odd photon to a level 
< 0.1%. 

To reject background from channels with muons or from possible misreconstruction 
of charged tracks, the muon detector is required to have no hits. 

Finally, the reconstructed kaon energy must be in the range 54 < Ek < 66 GeV 
and the reconstructed kaon mass be within 10 MeV/c^ of the PDG mass. 

With this selection, LI and L2 trigger efficiencies have been measured using control 
data samples recorded through minimum bias triggers. The LI trigger efficiency depen- 
dency on the cluster energies and relative positions has been studied using events with 
three photons in the calorimeter. As a result, a cut has been added in the selection of 
n^n'^'j events, which requires a minimum distance between photon clusters in both 
the X and y coordinates, as seen by the trigger projections. The LI trigger inefficiency 
has been measured as a function of the minimum photon energies. In order to exclude 
the region with smallest efficiency, a cut on the minimum photon energy of 5 GeV has 
been added to the event selection. After rejecting periods affected by identified hardware 
problems, the LI trigger efficiency is measured to be greater than 99%. 

The effect of the L2 trigger on the T* distribution has been studied using data and 
Monte Carlo simulation. The second level trigger applied an effective rejection of events 
with T* > 90 MeV, designed to provide high efficiency for the — > tc^tt^tc^ mode and to 
suppress — 7- TT^Tf^ decays. In order to take into account on-line resolution effects and 
avoid relevant biases, the corresponding cut has to be tightened in the off-line selection of 
7r^7r°7 to T* < 80 MeV. Data shows that the final L2 trigger efficiency is > 97% 
for all data taking periods, and is compatible with being flat in W. MC simulation of the 
L2 trigger efficiency shows that it is flat in W when a < 80 MeV cut is applied. Any 
departure from flatness is observed at low W only for events with 80 < < 90 MeV. 

4 Background contamination 

Due to their large branching fractions and their particle content, the most important 
potential background contributions to — )• t^^tt^'j come from — )• n^ir^n^ and 
— )■ 7r^7r° decays. In most of the previous experiments these backgrounds have been 
eliminated by demanding 55 < < 90 MeV ^\ As already mentioned, in NA48/2 a 
tighter cut of T* < 80 MeV has been implemented. Any remaining background coming 
from — )■ 7r^7r° would manifest itself as an excess of events in the region T*o ~ 110 MeV, 
T*o being the kinetic energy of the neutral pion in the kaon rest frame. No signiflcant excess 
of — tt^tt" events has been seen in the analyzed data sample, which is supported by 
the fact that data is compatible with signal in the high mass region (Figure [2]). 

The lower T* > 55 MeV cut is very efficient against -> Tr^Tr'^vr^ decays, but it 
also cuts away ~ 50% of the — )■ n^n^'j DE component. In this analysis an alternative 
procedure has been used in order to suppress the — )• Tr^Tr^vr'^ background. A — >■ 
7r^7r°7r° event can be reconstructed as — > 7r^7r°7, either if one of the four final state 
photons from the two vr^s is undetected (being emitted outside of the acceptance or being 
very soft) or if two of the photons overlap in the LKr calorimeter. These are mainly 
geometric effects and can be appropriately studied using Monte Carlo simulation. — )• 
n^n^n^ decays with one very soft photon in the final state are strongly suppressed by 

1) This cut is effective since T;{K^ n^n^) -110 MeV, while T;{K'^ -> tt^tt^tt^) <53 MeV. 



5 



phase space. 

In order to reject background from K"^ — )■ tt'^tx^tx^ decays with one undetected 
photon, two cuts have been implemented: the reconstructed kaon mass should be within 
±10 MeV/c^ of the nominal kaon mass and the center of energy (COG) of the system 
should be within 2 cm of the beam line. The center of energy is defined as the energy 
weighted radial position, at the LKr calorimeter front plane, of the decay particles in the 
event, the charged track being projected onto the LKr using the track directions before 
the magnet deflection. 

— )■ 7r^7r°7r° events with two photons overlapping in the LKr have the same 
signature as the signal in all detectors. Moreover, the reconstructed values of kaon mass 
and COG satisfy the selection criteria, as there is no energy loss for this topology. 

An algorithm has been developed to test the overlap hypothesis for all three photon 
clusters in the event. Let us consider the three reconstructed photon clusters with energies 
El, E2, and assume that the first one, with energy £"1, is really the overlap of two 
photons of energies E = xEi and E' = {1 — x)Ei. Assuming that the photon with energy 
E comes from the decay of the same vr*^ as the second photon cluster of energy E2, then 
the decay vertex for that 7r° would be: 

, _ v/(disti,2)^^^2 _ v/(disti,2)^a;Ei^2 , . 

where disti_2 is the radial distance between the photon clusters at the front of face of the 
LKr calorimeter, and M^o is the nominal 7r° mass. 

Similarly, assuming that the photon with energy E' comes from the decay of the same 7r° 
as the third cluster of energy E3, the decay vertex for this second vr" is given by: 

2 v/(disti,3)2E'£3 V(disti,3)'(l-a;)£i£3 



As the two neutral pions originate from the same kaon decay, they must satisfy the 

1 • 1 ^ 2 overlap 

constramt z = z 0= z n 

■j^U 

Hence Equations ([7]) and ([HD can be solved to obtain ^i^q"^^'"^. The procedure is 
repeated for all three photon clusters. The event is rejected if there is at least one solution 
in which ^"^'^'^''^^ jg compatible with zch within 400 cm. This procedure rejects ~ 98% of 
remammg — Tc'^Ti^n^ background to be compared with ~ 99% of T* >55 MeV cut 
and allows placing the lower T* cut at zero, increasing the sensitivity to the DE and INT 
components. 

Backgrounds from K^s and -R'^s, also accompanied by either a radiative or an ac- 
cidental photon, have been shown to be negligible using MC and the measured rate of 
accidental photons. 

The mass distribution of the K"^ — )■ '7r^7r°7 candidates is shown in Figure [2j All the 
background to K"^ — )• h'^ti^'j decays can be attributed to K"^ — )■ 7r^7r°7r° decays only, and 
is negligible in the selected mass region (< 10^^ with respect to the signal). 

5 Simulation: resolution and acceptances 

The NA48/2 Monte Carlo simulation is based on GEANT3 |23j. The K"^ 7r±7r°7 
matrix elements for IB, DE and INT [21 El [12] have been implemented separately, so that 
the three terms could be studied independently. No form factor for the DE contribution 
is considered at this stage. 
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Figure 2: Mass distribution of K"^ — )■ 71^71^7 reconstructed candidates. Simulated — )■ 
7r^7r°7 events are plotted in blue, while simulated — )■ vr^'^vr'^Tr'' events are superimposed 
in green. The data distribution (dots with error bars) is compatible with the sum of these 
two contributions, shown in red. Arrows delimit the selection region of |m^±^o^ — Mx| < 
10 MeV/(?, where the amount of background from K"^ — > vr^vr^vr^ decays is negligible 
with respect to the total sample. 

The complete detector geometry and material description is included in the sim- 
ulation. The energy deposited in the detectors during tracking of the decay particles is 
digitized and the reconstruction is performed as for real data. Detector imperfections, like 
DCH wire inefficiency maps, bad LKr calorimeter cells and malfunctioning readout cards 
are measured from data and included in the simulation in a run dependent way, so that 
experimental conditions are accurately reproduced. This procedure allows the L2 trigger 
to be reliably simulated. 

Non-gaussian tails on the cluster energy measurements may appear due to photon 
interaction with the LKr nuclei with a probability of 3x 10~^, resulting in an underestimate 
of the energy deposited in the electromagnetic shower. These tails have been parameterized 
from data events and a correction has been implemented for simulated cluster energies, 
modeling their effect. 

Using simulated events the relative resolution of W has been found to be ~ 1% in 
the range Q.2 <W < 0.9. 

To model the effect of the LI trigger, the measured efficiency shape as function of 
the minimum photon energy in the event has been implemented in the simulation. After 
requiring the photon energy to be greater than 5 GeV (see section |3]) the effect of the LI 
trigger was only significant for small W values, at the level of 1-2%. 

After all cuts the final acceptances for IB, DE and INT were respectively 3.13 %, 
4.45 % and 4.21 % in the range Q.2 <W < 0.9. 

6 DE and INT contributions 

The IB, DE and INT contributions have been extracted performing a fit of the data 
sample. Details on the fitting procedure, systematic uncertainties and final results are 
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given in the following. 

6.1 Fitting procedure 

A fitting program has been developed, based on a Poissonian Maximum Likelihood 
(ML) method, for the extraction of the IB, DE and INT contributions present in the data 
sample. The input information consists of the reconstructed W distributions of data and 
simulated IB, DE and INT samples. For each of them, 14 bins have been considered in 
the range 0.2 < W < 0.9. The program calculates the relative contribution of each com- 
ponent by minimizing the difference between the number of data events and the resulting 
total number of simulated events in every VT-bin. The fitted fractions are corrected for 
acceptance and final results are the relative contributions of DE and INT with respect to 
IB for T; < 80 MeV and < < 1: 

= (3.32 ± 0.15,^,4 ± 0.14,,,) X 10-2, (9) 

(-2.35 ± 0.35stat ± 0.39,,,) x lO^^. (10) 

The W distributions used in the fit contain 599k data events, and 3.770M simulated 
events (3.339M IB, 220k DE, 211k INT) in the fitting range 0.2 <W < 0.9. The correlation 
coefficient between the DE and INT fractions is —0.93. 

Figure E] a) shows the W data distribution superimposed with the simulated IB, 
DE and INT W distributions. The background contribution is also shown. The residuals 
of the ML fit to the data are shown in Figure [3] b). The of the residuals is 14.4 for 13 
degrees of freedom corresponding to a probability of Prob(x2)=0.35. 

Contour plots have been computed requiring the logarithm of the likelihood to 
change by 1.15, 3.1 and 5.9 units. These correspond to probabilities in the DE-INT frac- 
tions plane of 68.3 %, 95.5 % and 99.7 % respectively. They are shown in Figure HI 
including only statistical uncertainties. 

Figure [5] shows the ratio of the W distributions of data and IB simulation after all 
selection cuts and corrections. For large values of W the effect of the DE contribution is 
clearly seen. 

A fit to this ratio can be used to determine a and b as defined in Equation [3] and 
crosscheck the results obtained by the ML technique. Neglecting differences in acceptance 
distributions for IB, DE and INT as a function of W, evaluating the necessary integrals 
over the Dalitz plot, and fitting in the range 0.2< W <0.9, the result is: FracoE = 
(3.19 ± O.lQstat) X 10-2 and FracjNT = (-2.21 ± OAlstat) x 10"^ which agrees with the 
result obtained with the ML method (Equations [HI [ID]) within the systematic uncertainty. 

6.2 Systematic uncertainties 

The stability of the maximum likelihood fit result with respect to acceptance knowl- 
edge, misreconstruction effects and residual background contamination has been checked 
by varying the values of the main selection cuts within a reasonable range. The maximum 
likelihood fit has been repeated for every set of data and simulated samples obtained with 
the modified cut value, in order to study the differences in the resulting DE and INT 
fractions with respect to the standard (reference) ones. 
— Acceptance control. 

No significant effect was found varying the COG cut between 1 and 5.5 cm or the 
maximum distance between charged and neutral vertex \Az\ = \z.,^o — zcb\ between 
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(a) (b) 

Figure 3: (a) W distribution for data (squares) and sum of simulated signal and back- 
ground contributions (green). The lines show the signal contributions: IB (continuous), 
DE (dashed), and absolute value of INT term (dotted). Background contribution is also 
shown (red), (b) Residuals of the data with respect to the sum of simulated decay com- 
ponents weighted according to the Maximum Likelihood fit result. 



200 and 650 cm. The minimum distance between the center of the photon cluster 
and the pion impact point at the LKr was varied between 20 and 50 cm, and the 
results are in agreement within the uncorrelated uncertainties. 
The requirement that the photons cannot come from the vicinity of the DCHl inner 
flange was released and no change was observed in the result within the uncorrelated 
uncertainties. The acceptance definition of DCH 1, 2 and 4 was tightened by allowing 
the minimum and maximum radius of the track impact point to vary by a few cm 
away from standard cut values (12 cm and 150 cm). The results for DE and INT 
fractions did not change within the uncorrelated uncertainties. The distance between 
the neutral and charged vertex of the non-selected photon pairings was varied from 
to 800 cm, giving statistically compatible results. 

Limits on the contributions of the acceptance to the systematic uncertainties were 
estimated to be 0.10 x 10~^ for the DE term and 0.15 x 10~^ for the INT term. These 
values reflect the precision of the above studies involving the variation of the main 
selection criteria defining the acceptance. A confirmation of this estimate is given by 
the fact that the result of the a fit to the polynomial form (Equation [3] and Figure [5]), 
which uses the acceptance in a slightly different way, does not deviate from the main 
ML result by more than the previously quoted systematic uncertainties. 

— Background control. 

The E/p requirement for the charged track was varied between 0.75 and 1, with neg- 
ligible effect observed in the result. A 10 GeV/c cut is imposed on track momentum 
to maintain a high efficiency for the muon veto response, needed for background 
rejection. The value of this cut was varied between 5 and 15 GeV/c, giving statisti- 
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Figure 4: Contour plot for DE and INT terms. The cross shows the la statistical uncer- 
tainties on the projections. 

cally compatible results for DE and INT fractions. The reconstructed kaon mass is 
required to be within ±10 MeV/c^ of the nominal kaon mass, corresponding to ~ 4 
standard deviations a of the resolution. This value has been varied between 2.5 and 
9.5 a. It is expected that releasing this cut allows background from — ?> vr^vr'^Tr^ 
decays to enter the sample. However, the background is well enough separated from 
the signal, and no effect can be seen up to ~ 7 a. Therefore no systematic uncer- 
tainty was assigned due to background contamination. 
In addition to varying the relevant cuts, other effects have been studied. 

— Trigger efficiency. 

The effect of the LI trigger correction has been evaluated performing 1000 different 
fits. In each fit the bin values of the correction have been randomly varied according 
their uncertainties. Prom the distributions of fit results, rms values of 0.01 xl0~^ for 
DE and O.OSxlO'^ for INT have been extracted. 

The effect of a possible residual L2 trigger inefficiency has been studied changing 
the upper T* limit from 80 MeV to 65 MeV in five steps. To make a comparison with 
the standard measurement, all results have been extrapolated to the T* < 80 MeV 
region. A significant change in the result is observed only for the INT term. This 
behavior is in agreement with the hypothesis of a resolution effect of the L2 trigger 
response. The maximum observed difference, 0.3, has been assigned as systematic 
uncertainty to the INT result. 

— LKr energy reconstruction. 

The absolute LKr energy scale is known to a precision of 0.1 %. In order to evaluate 
the effect of this uncertainty, the data cluster energies have been multiplied by 
1.001. The results of the ML fit changed by 0.09x10"^ for the DE fraction and by 
0.21x10"^ for the fraction of INT. This difference has been assigned as systematic 
uncertainty to the result. 

Due to the zero-supression threshold applied to the LKr calorimeter cells at the 
readout, a non linear relation develops between the value of the energy deposited in 
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Figure 5: Ratio of W distributions for data events with respect to IB simulated events. 
The red hne shows the result of the fit to the polynomial form in Equation |3] in the range 
0.2< W <0.9. 

the calorimeter and its actual measurement. This non linear response is relatively 
more important for small cluster energies. Using — )■ 7r°7r° decays from the 2002 
run a correction has been applied to clusters with energies smaller than 11 GeV. 
This can be used for 2003 and 2004 data, where the readout thresholds were kept 
at the same value as in 2002. The shape of this non-linearity correction has been 
changed within a reasonable range. The results of the ML fit did not change with 
respect to the reference within the uncorrelated uncertainties. 

— Fitting method. 

Combinations of independent simulated samples of IB, DE and INT have been used 
as fake data, and fitted against standard simulated samples used in the data fit. In 
this way possible biases in the fitting technique have been excluded. 

— Radiative corrections. 

The simulation of the Inner Bremsstrahlung component has been interfaced with 
the PHOTOS package [21] thus generating one or more radiative photons. Only 1% 
of the selected events in this sample had more than one radiated photon. The ratio 
of the W distribution of these multi-photon events to the W distribution of the 
standard IB events was found to be a linear function of W. This model was used in 
a simple Monte Carlo program to evaluate its impact on the DE and INT fractions. 
This was found to be much smaller than 0.01 x 10^^ and therefore negligible. 

— Resolution control. 

Differences in resolution between data and simulation can produce different dis- 
tortions of the W shape, potentially biasing the DE and INT measurements. The 
comparison of kaon and neutral pion mass distributions in data and simulation 
shows relative resolution differences smaller than 2%. In addition, W distributions 
of simulated events have been smeared by different amounts and used as fake data 
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to be fitted against standard Monte Carlo samples with no extra smearing. These 
studies show that the systematic bias on the extraction of the DE and INT frac- 
tions is neghgible if the difference of W relative resolution between data and MC is 
smaller than 5%. 

The main contributions to the systematic uncertainty are summarized in Table [TJ 



Source of systematic uncertainty 


Effect on DE term 


Effect on INT term 


Detector acceptance 
LI trigger efficiency 
L2 trigger efficiency 
LKr Energy Scale 


< 0.10 X 10-^ 
0.01 X 10-2 
0.00 X 10-2 
0.09 X 10-2 


< 0.15 X 10-2 
0.03 X 10-2 
0.30 X 10-2 
0.21 X 10-2 


Total systematic uncertainty 


0.14 X 10-2 


0.39 X 10-2 



Table 1: Summary of non negligible systematic uncertainties. 



The effects of a possible form factor on the DE and INT terms measurements have 
been investigated. Following the work of reference [12] , a possible form factor in the pole 
part of the magnetic amplitude can be parameterized by the quantity rjv- The presence 
of such a form factor, if neglected in the analysis, induces anti-correlated variations of the 
extracted DE and INT terms such that AFracDE = — AFraciNT = — 0.01?7y Therefore the 
correction to the values of the DE and INT terms can be easily evaluated for any value 
of 7]v in the allowed range [0,1.5]. One should note that neglecting in the analysis the 
existence of a non- zero form factor present in the data, only a positive INT term could 
be induced, in no way faking a negative INT term. 

6.3 Results and discussion 

As seen from Equation [21 the simultaneous measurement of the DE and INT terms 
in K"^ — )■ 7r^7r*^7 decays allows to quantify the electric Xe and magnetic Xm contributions. 
Assuming a negligible amount of CP violation in the — )■ tt^tt^'j decays, the (p angle 
can be set to zero. The cosine of the difference between the two strong re-scattering phases 
can be approximated to 1 (see section [7]). The values obtained for Xe and Xm are 

Xe = (-24 ± Astat ± 4,,,) GeV-^ (11) 

Xm = (254 ± Gstat ± 6.,.) GeV'^ (12) 

The correlation coefficient is —0.83. 

The hypothesis that the chiral anomaly is the only source of magnetic amplitudes in 
the Direct Emission term predicts Xm ~ 270 GeV"^ [TJ [26] . Factorization models predict 
both an enhancement of Xm with respect to the pure chiral anomaly calculation and a 
positive Xe term in 7r±7r°7 [7l[Sl[27|. Our result shows a magnetic part compatible 

with the pure chiral anomaly and a negative Xe term, indicating that factorization models 
cannot provide an appropriate description of DE and INT in K"^ — t- tt^tc^'j. 

In order to compare the NA48/2 results with those from previous measurements, 
the ML fit of the selected sample has been redone setting the interference term to zero, 
and the result for DE extrapolated to 55 < < 90 MeV. 

Figure [6] shows the fit residuals. The demonstrates that the data distribution 
cannot be properly described without an interference term and that the DE-only fit is 
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Figure 6: Residuals of the Maximum Likelihood fit to the data with the INT contribution 
set to zero, with respect to the simulated distribution obtained from the weighted sum of 
the IB and DE components. 

not appropriate for this data. The result of this fit extrapolated to 55 < < 90 MeV is 
given here just for completeness: 

FracDE(INT = 0) = (0.89 ± QMstat ± 0.03,^,) x 10-^ (13) 

The systematic uncertainty assigned to this measurement corresponds to half of the dif- 
ference observed in the result varying the T* upper cut from 80 MeV to 65 MeV. 

Using the theoretical prediction for the branching ratio of the inner bremsstrahlung 
component as given in |3] and the current value for the branching ratio of the direct 
emission component as quoted in the PDG [2^ one obtains FracDE(INT = 0) = (4.3 ± 
0.7) X 10^7(2.61 X 10-^) = (1.65 ± 0.27) x lO^^. 

It should be noted that NA48/2 has been able to keep the rate of wrong solutions 
for the odd photon to < 0.1% while in other measurements [131 [IH lEl HSl |T7l [THl 
[inj the wrong solutions fraction was always larger than 10%. The relative background 
contamination is < 10~^ while in other experiments it is at the level of 10~^ or larger. 
The statistics used for this measurement is more than one order of magnitude larger than 
the sum of all previous experiments. In addition this measurement uses for the first time 
data without any lower cut on T*, with the benefit of better sensitivity to both DE and 
INT contributions to the K"^ — )■ 7r^7r°7 decay amplitude. 

7 CP violation 

As the decay of K"^ — )■ 7r^7r*^7 with direct photon emission is not suppressed by 
A/ = 1/2 rule, it has always been considered a good channel to search for CP violation 
[21 [2B]- According to |29j the magnitude of CP asymmetry in the Dalitz plot ranges from 
2 X 10~^ to 1 X 10~^ when the center of mass photon energy varies from 50 MeV to 
170 MeV. Possible supersymmetric contributions to direct CP violation in kaon decays 
can push the asymmetry to the level of 10"'^ in a specific region of the Dalitz plot |30j . 
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Present experimental knowledge on the asymmetry dates back to the late seventies and 
is (0.9 ±3.3)% [20]. 

In order to measure CP asymmetry in — )■ 7r^7r°7 decays, the previously described 
selection (Section [3]) has been modified to increase the statistical precision by setting the 
minimum photon energy to 3 GeV and eliminating the W > Q.2 cut. 

The new selection preserves the performance in terms of background rejection. Even 
if the selection criteria have not been completely optimized for studies of CP violating 
effects, it must be noted that the design of the experiment is such as to suppress beam 
and detector related differences between and K~ decays. 

For this reason differences for the two kaon charges in trigger efficiencies and ac- 
ceptances can be neglected, as proven in [21]. The residual effects are taken into account 
in the systematic uncertainty. To investigate CP violation in K"^ — 7r^7r°7, events were 
analyzed according to the reconstructed kaon charge, leading to a sample with 695k 
and 386k K~ . 

7.1 Integrated charge asymmetry 

The simplest observable that can be measured is the difference in the decay rates 
of and K~ . This can be expressed as the asymmetry on the total number of events 
Atv defined as: 

= itvm: 

where iV+, are the number of , K~ decays to ix'^-K^'-y in the data sample, and 
R is the ratio of the number of to K~ in the beam. Since R cannot be directly 
measured in NA48/2, normalization to another decay channel is required. For this work, 
the K"^ — )■ TT^TT^TT*^ decay has been chosen as normalization. The previously measured 
CP asymmetry in the rate for this channel is compatible with zero (-^ = (0.0 ± 0.6) 
10~^ PD]), and its uncertainty would contribute as an external source of uncertainty to 
our measurement of the CP rate asymmetry in the K"^ — )• 7r''=7r°7 channel. However, the 
NA48/2 data has shown no CP violation asymmetry in the Dalitz plot for K"^ — )■ 7r^7r''7r'' 
decays at the level of 10""^ [2T]. Since in most models the integrated rate asymmetries 
are expected to be smaller than the slope asymmetries (see for example [31]), we consider 
the absence of a CP rate asymmetry in the — t- vr^'^Tr'^Tr'' channel at the same level a 
plausible assumption, making the above external error negligible. 

A large number of — )■ 7r^7r''7r'' decays has been collected by NA48/2. Using the 
selection described in [52], only adapting few geometrical cut to fit the — )■ 7r''=7r°7 
ones, a measurement of the ratio R = 1.7998 ± 0.0004 has been performed with a high 
accuracy of 6R/R ~ 2 x 10""^. 

The ratio R has been computed in bins of kaon momentum and the corresponding 
asymmetry calculated. No dependency of the asymmetry was found neither as a function 
of T* nor as a function of the kaon momentum spectra. Differences in vr "'"-proton and 
TT^-proton cross section can induce a difference in the LI trigger efficiencies for and 
K~ events, that can refiect in a fake asymmetry in — )■ tt^tc^'j. The charged hodoscope 
inefficiency in detecting a pion depends on the pion hadronic cross section, and is different 
for the two charges. However, first order effects are also included in the R measurement 
performed using — )■ 7r''=7r''7r''. The major difference between — )■ ir^Tr^'y and — )■ 
7r^7r°7r° is due to different pion momentum spectra. The effect on the asymmetry has been 
calculated to be of the order of a few 10""^. A possible reconstruction induced asymmetry 
has been obtained to be < 5 x 10""^ using MC, while trigger effects have been evaluated 
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to be of the order of 4 x 10 ^. Finally, including the maximum R variation allowed by the 
NA48/2 estimation, the measured asymmetry is: 

Ar, = (0.0 ± 1.0,^,^ ± 0.6,,,) X lO-l (15) 

From the above value a limit for the rate asymmetry of I^tvI < 1.5 x 10"'^ at 90% 
CL can be deduced. An alternative approach would have been to normalize the K'^ and 
K~ fluxes using inclusive decay channels, such as the inclusive three pion channels or the 
K"^ — )■ vr^vr^ plus K"^ — 7r^7r°7 channels. For those the absence of a CP asymmetry is 
guaranteed by the CPT theorem together with unitarity when small final state interactions 
are neglected [MllS^. However, in the NA48/2 experiment other possible normalization 
channels different from — )■ vr^Tr^Tr" have been collected with a trigger different from 
that used for K"^ — )■ n'^'K^'y decays, so that using them in the asymmetry measurement 
would require an accurate knowledge of the relative trigger efficiencies. Therefore such an 
alternative is not presented at this time. 

The asymmetry A^- can be related to the CP violating phase (p by Equation [2j 

r+ — F+ — F^ 
An = p , , -p_ ~ = {Iint/ IiB)2XEmj^ml sin (psm{5l - 5^) = 

= e(W/^,B), (16) 



{IiNT/IiB)2XEmlmlsm{6l - 



where the parameter e, describing the asymmetry, has been introduced. 

The difference between the two strong re-scattering phases 6q is evaluated to be 
6.6° ± 0.5°. This is calculated using the theoretical predictions for 6l and 6^ as a function 
of the TTTT mass from [25], weighted by the observed Mtt+tto data distribution. 

Using in addition the value of 0.105 for the ratio of W integrals, and Xe from 
the NA48/2 measurement (Equation [TT]) a measurement for sin0 = —0.01 ± 0.43, or 
I sin 01 < 0.56 at 90% CL is obtained. 

7.2 Fit to the W spectrum 

Another interesting check of CP violation can be obtained by looking at the distribu- 
tion of the asymmetry as function of the Dalitz plot variable W. In fact, an enhancement 
of the asymmetry in particular regions of the Dalitz plot is suggested by [221 SO] • Using 
Equation [21 the following W dependence of the asymmetry (Aw) is predicted: 

W - ^(l + ("±^)«'= + (IB) 

dr+ dr- 

r 1/1/ - 

(19) 



dAw _ iw-lw _ eVT 



dW ^ + ^ l + aW^ + hW^' 

dW dW 

where a and h have been already extracted in section [6l neglecting CP violation effects. 
The parameter e is the only free parameter left in the fit. Figure [7] shows the asymmetry 
as a function of the W variable for the whole data sample together with the fit result. 
Multiplying the e parameter obtained from the fit with the value of the integral ratio 
IintIUb = 0.105, the value of Aw, the asymmetry in the W spectra, is: 

Aw = cIint/Iib = (-0.6 ± 1.0,^^0 x 10"^ (20) 

The value of Aw is compatible with the result An of the overall charge asymmetry. 
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Figure 7: Measured asymmetry as defined in Eq. [12] as a function of W. 



8 Summary 

From a sample of about 600k — )■ 7r^7r°7 decay candidates, the NA48/2 exper- 
iment has measured the relative amounts of DE and INT with respect to the internal 
bremsstrahlung (IB) contribution in this decay in the range < < 80 MeV: 

FracDE(0 < T; < 80 MeV) = (3.32 ± 0.15,tat ± 0.14,^,,) x 10-^, 
FraciNT(0 < T; < 80 MeV) = (-2.35 ± 0.35,tat ± 0.39,^,,) x lO^^. 

The relative background contamination has been kept to < 10~^, and the rate of wrong 
solutions for the odd photon to < 0.1%. Thanks to the implementation of an algorithm 
rejecting background from — )■ n'^n^n'^ decays, the cut on T* could be released below 
the standard 55 MeV used by most of the previous experiments, gaining in sensitivity to 
both DE and INT contributions to the K"^ — t- tt'^tt^'j decay amplitude. This measurement 
constitutes the first observation of an interference term in — )■ tt^tt^'j decays. The 
results for electric and magnetic contributions, Xe = (—24 ± Agtat ± 4:sys) GeV"'^ and 
Xm = (254 ± 6stat ± Qsys) GeV~^ indicate that factorization models cannot be applied to 
estimate the size of DE and INT in — )■ n'^TT^'y and that the magnetic part is compatible 
with pure chiral anomaly. 

Using a slightly modified event selection, two samples of 695k and 386k K~ 
have been reconstructed and used to set a limit on the CP violating asymmetry in the 
and K~ branching ratios for this channel of less than 1.5 xlO~^ at 90% confidence 
level. For this measurement the K"^ — )■ tt^tt^tt^ decay has been used for normalization. 
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